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_i 1 

The specific heat of liquid helium was calculated theoretically in the Landau 
theory ^. The results deviate from experimental data in the temperature 
region of 1.3 - 2.1 K. Many theorists subsequently improved the results of 
the Landau theory by applying temperature dependence of the elementary 

'^ ' excitation energy^,^. As well known, many-body system has a total energy 

of Galilean covariant form. Therefore, the total energy of liquid helium has a 
Q \ nonlinear form for the number distribution function. The function form can 

be determined using the excitation energy at zero temperature and the latent 
heat per helium atom at zero temperature. The nonlinear form produces new 
temperature dependence for the excitation energy from Bose condensate. We 

\^ \ evaluate the specific heat using iteration method. The calculation results of 

cn ■ the second iteration show good agreement with the experimental data in the 

temperature region of - 2.1 K, where we have only used the elementary 

y^ \ excitation energy at 1.1 K. 

OO 

O 



1. Nonlinear form of total energy 

Liquid iieliuiii system lias a total Hamiltonian as 

p2 1 

^ = Y1 ^"^P^P + 2F ^ 5(k)a;+kOq-kapaq (1) 

P p,q,k 

where m is the mass of a helium atom, a* and a^ respectively signify the cre- 
ation and annihilation operators. We examine the general form of the total 
energy via the unitary transformation U diagonalizing the Hamiltonian H. 
All eigenstates are described as | eigenstate >= Ua*a*a*---a*j^ \ > 
where [ > denotes the vacuum state. New creation and annihilation op- 
erators are defined as 



Ap = Ua*pU-^,Ap = UapU'^ (2) 

which indicate the creation and annihilation operators of a quasi-particle. 
We designate this quasi-particle as a "dressed boson". The dressed boson 
number operator is defined as 

np = A;Ap. (3) 

The total number conservation and the total momentum conservation are 
expressed as 



q q q q q 

That is to say, the total number of helium atoms is equal to the total number 
of dressed bosons and the total momentum of helium atoms is equal to the 
total momentum of dressed bosons. 

The total energy of the system is a sum of the kinetic energy K of the 

q2 

center of mass and Galilean invariant terms X : K = ^^ where M is the 
total mass of liquid helium. 

p p.q 

^2 



yj rip + (Galilean invariant terms) (6) 



P 2m 

where Galilean invariant terms are described only by relative momenta of 
dressed bosons: 



(Galilean invariant terms) = 
-TT Yl /2(P - q)npnq + ^ X! /3(P - q, P - k)npnqnk + . . . (7) 

P,q p,q,k 



Substitution of Eq.© into Eq.® yields 



p2 1 

^ = Y1 ^rr^p + 77 5Z /(p - q)'^p"q (8) 



P 2"^ ^ P,q 



where we abbreviate higher terms because three-particle colhsion is a rare 
case for diluteness of Hquid hehum. The single excitation state has a dis- 
tribution of {no = A^ — l,np = 1} and therefore its total energy is derived 
from Eq.Q as follows: 

^ = /(0)iV + |^ + 2(/(p)-/(0)) (9) 

where we have used 1/A^ ~ and the spherical symmetric property of the 
function /(p). Therein, the latent heat at zero Kelvin is equal to —Nf{0). 
Accordingly the elementary excitation energy at zero Kelvin is given by 

2 

^p~2m^ 2(/(p) — /(O)). This relation engenders a function form of the 
nonlinear term as-^ 



2. Coupled equation determining distribution of dressed bosons 

The energy of one dressed boson is an increase value of the total energy 
when one dressed boson is added to the system. Accordingly the dressed 
boson energy is defined as Wp = 6E/5np. The calculation result for the 
derivative of Eq. ([8]) shows 

^p^ = 1^ + 1 E /(p - q)"q - ]^ E /(« - t)^snt (11) 

q s,t 

where we have used /(p — q) = /(q — p)- The distribution function is 

determined as 

_ 1 

""P " exp((^p(r)-^)/(A;BT))-l ^^^^ 

We can obtain approximate solutions of the coupled equations of (fTTI) 
and ([12]) via the iteration method ™ We adopt the Landau distribution 
function as the zero-th order distribution: 

no{p,T) = (13) 

exp(e^/(A;BT)) - 1 

The j-th order solutions are derived from the (j-l)-th distribution function 
as follows: 

'^j(P'^) = |;:+^E/(p-qK-i(q'^)-42 E/(«-tK-i(«'^K-i(t'^) 



2m' N^-""^ Mr-j-ivM,^. ^2 

q s,t 



(14) 



This j-th energy form produces the j-th distribution function: 

1 



nj{p,T) 



exp((u;,(p,r) - io,{0,T))/ikBT)) - 1 



(15) 



Therein the excitation energy from the Bose-Einstein condensate of dressed 
bosons is expressed as 



(^j{p,T) =ujj{p,T) -ujj{0,T). 



(16) 



We can evaluate the second order solutions £2 {p, T) and n2 (j>, T) via the 
iteration processes from the zero-th order distribution. 



3. Evaluation of specific heat 

Using the second order excitation energy €2 {p, T) and the distribution 
function n2{p,T), we can calculate the second order approximation values 
of specific heat as follows ^, 



CP = TkBj^{^}^Io^ {log(l + n2(p,r)) + ^-2gln2{p,T)}pMp 



+ 



+ TtwIo°'in2ip,T)feMe2(j^,T)/{kBT)) 






e2(p,T) r de2{p,T) } I 2^^ 



ksT \ dT jp 

The evaluated results are shown in Fig[T] and Figl2l 
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Figure [H and [2] indicate the second order results of specific heat via the 
nonlinear theory. The curves express the calculated values. The dots with 



Cg [J/(mol«K)] 
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Fig. 3. The large dots colored with blue indicate the experimental data 
'^, and the small dots with red are measured by Lipa et al - . BD theory 
indicates the results in the reference oi'^. BCY theory indicates the results 
in the reference oi'^. The curve of the nonlinear theory is the results of the 
reference '^ 



red indicate experimental datal^. As shown in Fig(T]and Figl2l the theoretical 
values of the second order are in good agreement with the experimental 
data for T < 2. IK. It is noteworthy that the present calculation uses the 
experimental values of excitation energy only for the temperature 1.1 K. Of 
course the iteration method is insufficient in close vicinity of the A transition 
temperature. We have discussed origin of the logarithmic divergence at the 
A point in the previous paper'^. It is clarified that the logarithmic divergence 
is caused by the nonlinear mechanism of the total energy. The calculation 
results are shown in Fig. [3l 

Thus the nonlinear theory has well explained the temperature depen- 
dence of the specific heat of superfluid helium for all temperature region. 
Accordingly the nonlinear mechanism of total energy is important for un- 
derstanding the properties of liquid helium. 
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